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Abstract 

Mosquito-borne diseases present some of the greatest health challenges faced by the world today. In many cases, existing 
control measures are compromised by insecticide resistance, pathogen tolerance to drugs and the lack of effective vaccines. 
In light of these difficulties, new genetic tools for disease control programmes, based on the deployment of genetically 
modified mosquitoes, are seen as having great promise. Transgenic strains may be used to control disease transmission 
either by suppressing vector populations or by replacing susceptible with refractory genotypes. In practice, the fitness of 
the transgenic strain relative to natural mosquitoes will be a critical determinant of success. We previously described a 
transgenic strain of Anopheles gambiae expressing the Vida3 peptide into the female midgut following a blood-meal, which 
exhibited significant protection against malaria parasites. Here, we investigated the fitness of this strain relative to non- 
transgenic controls through comparisons of various life history traits. Experiments were designed, as far as possible, to 
equalize genetic backgrounds and heterogeneity such that fitness comparisons focussed on the presence and expression of 
the transgene cassette. We also employed reciprocal crosses to identify any fitness disturbance associated with inheritance 
of the transgene from either the male or female parent. We found no evidence that the presence or expression of the 
effector transgene or associated fluorescence markers caused any significant fitness cost in relation to larval mortality, pupal 
sex ratio, fecundity, hatch rate or longevity of blood-fed females. In fact, fecundity was increased in transgenic strains. We 
did, however, observe some fitness disturbances associated with the route of inheritance of the transgene. Maternal 
inheritance delayed male pupation whilst paternal inheritance increased adult longevity for both males and unfed females. 
Overall, in comparison to controls, there was no evidence of significant fitness costs associated with the presence or 
expression of transgenes in this strain. 
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Introduction 

Mosquito species are important global vectors of diseases such 
as malaria, dengue, yellow fever, filariasis and encephalitis. Of 
these, malaria is arguably the greatest problem, with up to 40% of 
the human population at risk of infection and around one million 
deaths each year. One of the major routes to disease control is to 
tackle the mosquito vectors themselves, primarily with insecticides. 
However, as a result of these interventions, mosquitoes are 
becoming increasingly resistant to pesticides [1] and additional 
approaches are clearly required. The feasibility of utilising 
transgenic mosquito technology to control disease transmission 
has been demonstrated [2-11]. However, as with any genetic 
intervention, the ultimate goal of disease control through 
transgenic technology will rely upon effective inheritance of the 
transgene and its introgression into field populations. For such 
technology to become a viable part of an integrated pest 
management programme, genetically modified mosquitoes must 
be able to compete successfully for available resources. This 
requires either the absence of a fitness load associated with the 
genetic modification, or else a modest fitness load that may require 



the use of a gene drive mechanism to help fix the transgene in the 
population [12,13]. To achieve this outcome it is important that 
transgenic mosquitoes, destined for mass rearing and field release, 
have minimal fitness costs associated with transformation and 
transgene expression. 

The overall fitness of a given transgenic strain will depend upon 
a number of factors that each require investigation prior to 
consideration for use in field trials. First, the introduction of a 
designed intervention against the pathogen may have unintended 
fitness consequences for transgenic strains. For example, the 
introduction or over-expression of effector genes may have a 
negative impact on fitness [14] and the expression of different 
effector genes (even under control of the same regulatory 
sequences) can have differing fitness costs in transgenic strains 
[15]. Secondly, factors that are independent of the effector gene or 
refractoriness to disease transmission may impact fitness. For 
example, fluorescence marker expression cassettes are typically 
employed to allow identification of transformed individuals [16]. 
There are reports of adverse fitness effects on transgenic 
mosquitoes caused by simple promoter-marker combinations 
[17,18]. This may depend upon the regulatory sequences chosen 
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for fluorophore expression but is most likely a result of insertional 
mutagenesis following transgene insertion. Insect transgenesis is 
typically mediated by transposons that have an essentially random 
pattern of genome integration [19]. There is therefore an 
associated potential to disrupt native gene function [14]. This 
can impose fitness costs, but is largely dependent on the choice of 
transposon. For example, the P element exhibits a strong tendency 
to integrate into gene promoter regions whereas the pigsyBac 
element appears to have evolved to largely avoid host gene 
disruption [20]. In addition, the essentially random nature of 
transposon integration can result in position effects whereby 
transgene(s) are influenced by the proximity of other gene 
regulatory sequences, or epigenetic modifications. In our experi- 
ence, position effects can have a significant impact on the observed 
levels of both fluorescence marker and effector gene expression 
[9,16,21]. It follows that transgene location is an important factor 
to consider and that fitness implications associated with transgene 
expression could be either ameliorated or compounded by 
genomic location. Thirdly, and perhaps most importantly, 
inbreeding and loss of genetic heterogeneity as a result of the 
genetic botdenecks inherent in genetic modification may impact 
on fitness. Transgenic strains are established from very few 
founder individuals (typically from a single transgenic adult) and 
are then maintained as small laboratory populations. As a result, 
the resulting strains may carry fixed deleterious alleles with the 
potential to impact on fitness [18]. Indeed, the detrimental effects 
of inbreeding on mosquitoes in the absence of transgenes has been 
previously demonstrated [13]. 

A range of fitness studies have been undertaken to determine 
the impact of transgenesis on mosquito biology in both Anopheles 
stephensi [6,10,12,15,18,22-24] and Aedes aegypti [13,25,26]. Signif- 
icant fitness disturbances have been reported for some transgenic 
strains but the effect is not uniform across similar strains and can 
be impacted by any or all of the factors discussed above. The 
effector gene itself has been directly responsible for an impact on 
fitness in two studies [6,15]. Corby-Harris et al. reduced malaria 
parasite infection by increasing Akt signalling in the mosquito 
midgut [6] . Lifespan was also reduced in these transgenics, even in 
non-blood-fed females, which was attributed to leaky expression 
from the carboxypeptidase promoter. A reduction in fitness was 
also observed when phospholipase A2 from honeybee venom was 
expressed from the carboxypeptidase promoter [15]. This was 
presumed to be as a consequence of internal damage to the midgut 
and was resolved by expression of a mutant protein [7] . The same 
researchers observed no fitness disturbance in hemizygous 
transgenic mosquitoes expressing an alternative antimalarial 
effector gene, SMI, in the midgut [15]. A fitness advantage 
compared to wild-type was even reported when SMI transgenic 
mosquitoes were maintained on Plasmodium infected blood [22]. 
Isaacs et al. [10] also reported no fitness load associated with the 
expression of single-chain antibodies when compared to the 
control strain. Since females of one strain exhibited an increased 
median lifespan, there may once again be an associated fitness 
advantage. Where a fitness load has been identified, it has more 
often been attributed to position effects or insertional mutagenesis, 
compounded by inbreeding, rather than expression of the 
transgene itself. Such a fitness load was reported for three 
independent, homozygous transgenic strains expressing SMI into 
the haemocoel [12]. Due to the particular inbreeding issues 
associated with generating homozygous transgenic strains, which 
are likely to exacerbate founder effects and result in inbreeding 
depression, it is not surprising that fitness loads have also been 
identified in other homozygous strains [18,25,26]. Such fitness 
disturbance is more pronounced when assessing changes in 



transgenic allele frequencies over time rather than measuring 
fitness parameters directly [12]. Fitness loads associated with 
position effects are dependent upon transgene location and result 
in a lack of uniformity between strains such as those reported for 
fitness profiles of four independent hemizygous An. stephensi 
docking strains [23]. However, this study concluded that none of 
the strains carried significant fitness burdens compared to their 
non-transgenic counterparts. Other studies, designed to separate 
the effect of inbreeding from that of genetic modification have 
concluded that transgenesis has an effect on fitness that is 
additional to that of inbreeding. Koenraadt et al. [13] found that 
transgenic Ae. aegypti fared least well in comparison to both wild- 
type and an inbred strain for most fitness parameters measured, 
but that the inbred strain was also less fit than wild-type. Another 
study in which transgenic Ae. aegypti were outcrossed to the 
unmodified strain for five generations in order to minimise 
inbreeding effects, observed significant differences in all measured 
fitness parameters when comparing the outcrossed transgenic 
strain to wild- type [13]. It was therefore suggested that in this 
instance transgene products may be responsible for the overall 
reduction in fitness of the transgenic strain compared to wild-type. 

It is clear from these various fitness studies that, in designing 
comparisons between transgenic and non-transgenic strains, care 
needs to be taken to equalize as far as possible the genetic 
backgrounds and heterogeneity of the strains under comparison. 
Only then can any inference be drawn about the possible impact 
of the genetic modification itself. The present study was designed 
to investigate fundamental fitness parameters of the transgenic An. 
gambiae EVida3 strain [9]. This strain expresses a synthetic 
antimicrobial peptide, Vida3 [27], into the female midgut lumen 
following a blood-meal. EVida3 was generated using the phiC31 
site-specific integration system [28], whereby the effector gene 
cassette was integrated into a previously characterised docking site 
on chromosome 3R. Therefore, in addition to the Vida3 cassette, 
the strain expresses two fluorescence markers (ECFP and DsRed2), 
both under control of the eye-specific 3 xP3 promoter. The fitness 
costs associated with inbreeding have been well documented 
[12,13] and a study modelling the population dynamic effect 
following the introduction of a genetically modified mosquito into 
a wild-type population concluded that the relative fitness of the 
resulting hemizygotes is crucial for a successful release [29]. 
Hemizygotes were therefore chosen for this study and were 
generated by outcrossing EVida3 homozygotes to the parental 
KIL laboratory strain from which they were derived in order to 
equalize genetic backgrounds [25]. Hemizygotes were established 
using reciprocal crosses between EVida3 and KIL to enable 
comparisons between strains that inherited the transgene via the 
maternal parent (EV3M) or paternal parent (EV3P). This 
assessment of fitness parameters is an important step towards 
increasing our understanding of the biology of transgenic strains 
such as EVida3 and their potential performance in the field. 

Results 

Larval Mortality 

Control (KIL) and hemizygous transgenic (EV3P and EV3M) 
larvae were generated and reared following detailed protocols to 
minimise the impact of environmental variance. Three indepen- 
dent experiments were established, each comprising 15 replicates 
with 30 individual mosquitoes per strain. Replicates were 
individually randomised within the insectary to minimise con- 
founding placement effects. Larval death was recorded at the same 
time each day until all larvae had either pupated or died. Mean 
larval mortalities were 2.69 (9.0%), 0.96 (3.2%) and 1.67 (5.6%) 
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Figure 1. Larval mortality, pupal sex ratios and wing lengths. Data represent comparisons between the control (KIL) and hemizygous 
transgenic strains (EV3P and EV3M). (A) Larval mortality data are pooled over 3 independent experiments, each with 15 replicates of 30 individuals per 
strain. Boxplots show the interquartile range (open boxes), the median (horizontal line for KIL but coinciding with the interquartile range at 1 for both 
EV3P and EV3M) and the maximum and minimum values (vertical lines). There were no significant differences between strains (p = 0.368). (B) 
Numbers of pupae and pupal sex ratio data are pooled over 3 independent experiments, each with 15 replicates of 30 individuals per strain. The 
percentage of females (open bars) and males (shaded bars) for each strain is indicated. Statistical analysis showed no significant difference in sex 
ratios between strains (p = 0.673). (C) Wing length data are pooled over 3 independent experiments, each involving 100 female and 100 male 
mosquitoes for each strain. (Ci) Female wing lengths were non-normally distributed and are represented as boxplots showing the median (horizontal 
line), interquartile range (boxed) and maximum and minimum values (vertical lines). There were no significant differences between strains (p = 0.1 12). 
(Cii) Male wing lengths followed a normal distribution and are represented as open bars with standard errors (vertical lines). There were no significant 
differences between strains (p = 0.590). 
doi:10.1371/journal.pone.0088625.g001 



for KIL, EV3P and EV3M respectively. Since these data were 
non-normally distributed, the medians with associated ranges were 
plotted (Figure 1A). Statistical analysis using a non-parametric 
Kruskal-Wallis test revealed no significant difference in larval 
mortality between strains (p = 0.368). 

Pupal Sex Ratio 

The number and sex of all pupae recovered from the larval 
mortality experiment described above were analysed using data 
pooled across all replicates from 3 independent experiments. From 
a total of 1350 larvae of each strain, 1229 (91%), 1307 (96.8%) and 
1275 (94.4%) larvae survived to pupation for KIL, EV3P and 
EV3M, respectively. The sex ratios (female:male) were 48.7:51.3 
(KIL), 50.5:49.5 (EV3P) and 49.8:50.2 (EV3M), respectively 
(Figure IB). These data were analysed using a Chi-squared 



contingency test which revealed no significant difference in sex 
ratio between the three strains (% 2 = 0.792; 2 degrees of freedom; 
p = 0.673). 

Wing Length 

Mosquitoes for wing length measurements were collected from 
the longevity study described below. Wing lengths were measured 
for both males and females from control (KIL) and hemizygous 
transgenic (EV3P and EV3M) strains. Data were pooled across 3 
independent experiments each comprising 100 males and 100 
females of each strain. Average female wing lengths for each strain 
were 3.11 mm (KIL), 3.16 mm (EV3P) and 3.16 mm (EV3M) 
respectively (Figure 1 Ci). Female wing length data were non- 
normally distributed (DAgostino-Pearson omnibus test) and were 
compared using the Kruskal-Wallis test, which revealed no 
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significant difference between strains (p = 0.112). Average male 
wing lengths for each strain were 2.95 mm (KIL), 2.96 mm 
(EV3P) and 2.95 mm (EV3M) respectively (Figure ICii). Male 
wing lengths followed a normal distribution (D'Agostino-Pearson 
omnibus test) and were analysed using ANOVA which revealed no 
significant difference between strains (p = 0.590). 

Age at Pupation 

Data for age at pupation for control (KIL) and hemizygous 
transgenic (EV3P and EV3M) strains were collected from the 
larval mortality experiment described above. Age at pupation was 
recorded separately for males and females of each mosquito strain. 
For females, the mean age at pupation was 7.76 days (KIL), 7.72 
days (EV3P) and 7.76 days (EV3M) respectively. Statistical 
analysis of pooled data using the non-parametric Kruskal-Wallis 
test revealed no significant difference in age at pupation between 
strains (p = 0.367; Figure 2A). For males, the mean age at pupation 
was 7.66 days (KIL), 7.59 days (EV3P) and 8.32 days (EV3M) 
respectively. Statistical analysis of pooled data (Kruskal-Wallis test) 
revealed a significant difference in age at pupation between strains 
(/><0.0001; Figure 2B). Further analysis using Dunn's multiple 
post-test for pair-wise comparisons identified a significant differ- 
ence between the age at pupation of E V3M males when compared 
to either KIL or EV3P (p<0.0001 by Mann Whitney test in both 
cases; Figure 2B). There was no significant difference in age at 
pupation between KIL and EV3P males (Figure 2B). 

Adult Longevity 

Mosquitoes were reared following detailed protocols to mini- 
mise environmental variance and adult longevity was recorded 
separately for males and females of each strain. Data were pooled 
from 3 independent experiments, each comprising 100 male and 
100 female mosquitoes per strain, separated into 10 replicates of 
10 newly eclosed mosquitoes. Adults were reared on sugar and 
water and dead mosquitoes were removed and recorded daily until 



all mosquitoes had died. For females, average adult longevity was 
28.57 days (KIL), 34.32 days (EV3P) and 30.37 days (EV3M) 
respectively. For males, average adult longevity was 29.50 days 
(KIL), 33.16 days (EV3P) and 26.24 days (EV3M) respectively. 
Pooled data were analysed for statistical significance using the log- 
rank (Mantel-Cox) test. For both KIL and EV3P, adult longevity 
did not differ significantly between males and females (p = 0.1649 
and p = 0.0585 respectively; compare Figures 3A, 3E) whereas for 
EV3M, females were significantly longer lived than males 
{p — 0.0023; compare Figures 3A, 3E). Further analysis revealed 
a significant difference in longevity between strains for both 
females (Figure 3A; /K0.0001) and males (Figure 3E; /K0.0001). 
Pooled female and male data were then separately subjected to 
pair-wise comparisons to determine where the major differences in 
longevity lay. For females, there was no significant difference in 
longevity between KIL and EV3M (Figure 3C; ^ = 0.462) but 
significant differences were observed between KIL and EV3P 
(Figure 3B; /K0.0001) and between EV3P and EV3M (Figure 3D; 
/><0.0001). For males, significant differences were identified in all 
pair-wise comparisons, between KIL and EV3P (Figure 3F; 
p = 0.004), between KIL and EV3M (Figure 3G; p = 0.001) and 
between EV3P and EV3M (Figure 3H; /K0.0001). 

Fecundity and Hatch Rate 

Mosquitoes for fecundity and hatch rate assessments were 
reared together with those for the blood-feeding longevity 
experiments detailed below, following detailed protocols to 
minimise environmental variance. Representative samples from 
each of 3 independent experiments were measured for wing 
length. Average female wing lengths for each strain were 3.29 mm 
(KIL), 3.30 mm (EV3P) and 3.34 mm (EV3M) respectively 
(Figure 4A), which were larger than those reared previously 
(Figure ICi). These non-normally distributed data (DAgostino- 
Pearson omnibus test) identified a significant difference in wing 
length (p = 0.0013) using the Kruskal-Wallis test for this set of 
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Figure 2. Age at pupation. Data represent comparisons between the control (KIL) and hemizygous transgenic strains (EV3P and EV3M) and show 
age at pupation for females (A) and males (B). Data are pooled over 3 independent experiments, each with 1 5 replicates of 30 individuals per strain. 
Boxplots show the interquartile range (open boxes) with maximum and minimum values (vertical lines). Different lower case letters (a, b) identify 
significant differences. The median pupation age is 8 days for all data sets. The age at pupation for females was not significantly different between 
strains (p = 0.367) but there was a significant difference in male pupation rates between strains (p<0.0001). Further analysis identified a significant 
difference between EV3M and KIL (p<0.0001) and between EV3M and EV3P (p<0.0001). 
doi:1 0.1 371 /journal.pone.0088625.g002 
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Figure 3. Adult longevity survival curves. Graphs represent comparisons between the control (KIL, black) and hemizygous transgenic strains 
(EV3P, green and EV3M, red) and show Kaplan-Meier adult survival curves for females (A-D) and males (E-H). (A) shows survival curves for females of 
all three strains. (B-D) show paired comparisons between strains for female survival. (E) shows survival curves for males of all three strains. (F-H) show 
paired comparisons between strains for male survival. All survival curves represent pooled data from three independent experiments, each with 10 
replicates of 10 adults for each strain. Significance levels are displayed on each graph (***, p<0.0001; * p<0.05; ns; not significant). 
doi:1 0.1 371 /journal.pone.0088625.g003 



samples. Pair- wise analysis (Mann- Whitney test) showed that 
EV3M females were significantly larger than both EV3P 
(p = 0.01 1) and KIL (p = 0.0003), but that there was no significant 
difference in size between EV3P and KIL (p = 0.507). Egg batch 
sizes collected from individual females followed a normal 



distribution (D'Agostino-Pearson omnibus test) and were analysed 
using ANOVA, which also revealed significant differences between 
strains (/)<0.0001; Figure 4B). Mean egg batch sizes were 129.1 
(KIL), 142.7 (EV3P) and 164.2 (EV3M) respectively and pair-wise 
analysis (unpaired t-test with Welch's correction) identified 
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significant differences between all 3 strains. EV3M egg batches 
were significantly larger than both EV3P (0 = 0.0025) and KIL 
(0<O.OOO1), which is most likely attributable to the larger wing 
length of EV3M females (Figure 4A). A significant difference was 
also identified between KIL and EV3P egg batches (0 = 0.0411), 
where EV3P produced slighdy larger egg batches than the control, 
although their wing lengths were not significantly different 
(Figure 4A). The mean hatch rates for egg batches collected 
above were 85.72% (KIL), 86.79% (EV3P) and 86.36% (EV3M; 
Figure 4C). Statistical analysis of pooled data (Kruskal-Wallis test) 
revealed no significant differences in hatch rates between strains 
(0 = 0.773). 

Female Longevity Following Blood-feeding 

Control and hemizygous transgenic mosquitoes were reared as 
detailed in Materials and Methods and adults collected daily. 
Three days post-eclosion, female sub-populations were either left 
unfed or offered a blood-meal (designed to stimulate expression of 
the Vida3 peptide from the carboxypeptidase promoter), from 
which fully engorged females were selected. Blood-fed cages were 
offered oviposition sites two and three days post-blood-meal and a 
second blood-meal eight days post-eclosion. Dead mosquitoes 
were removed and recorded daily and data were pooled from 3 
independent experiments, each comprising 4x25 non-fed and 
4x25 blood-fed mosquitoes per strain. For mated but non-blood- 
fed females average longevity was 12.90 days (KIL), 14.47 days 
(EV3P) and 12.55 days (EV3M), which is considerably shorter 
than for females reared previously (Figure 3), which were neither 
mated nor subjected to periods of sugar starvation. Data were 
analysed for statistical significance using the log-rank (Mantel-Cox) 
test which revealed a significant difference between strains 
(0 = 0.0004; Figure 5A). Separate pair-wise comparisons identified 
unfed EV3P females to be significandy longer lived than both KIL 
(0 = 0.002) and EV3M (0 = 0.0002), although there was no 
significant difference in longevity between unfed KIL and 
EV3M (0 = 0.524). For females that were blood-fed, more than 
95% of surviving individuals from all three strains took a second 
blood-meal five days later. Average longevity for blood-fed females 
was 15.81 days (KIL), 16.07 days (EV3P) and 16.68 days (EV3M). 
Analysis by log-rank (Mantel-Cox) revealed no significant differ- 
ence between blood-fed females (0 = 0.1 54; Figure 5B). Comparing 
Figures 5A and 5B, there was a significant increase in longevity in 
sugar-starved females that were blood-fed compared to non-fed for 
both KIL and EV3M (0<O.OOO1 for both) but not EV3P 
(0 = 0.159). 

Discussion 

This study compared the fitness of hemizygous transgenic 
mosquitoes expressing an antimalarial effector gene [9] to the 
wild-type host strain from which they were derived. The 
hemizygous strains were established from reciprocal crosses 
between KIL and EVida3 homozygotes such that the transgene 
was inherited from either the paternal (EV3P) or the maternal 
(EV3M) parent. Hemizygotes were used in an attempt to equalize 
the genetic backgrounds and heterogeneity of the strains under 
comparison and to minimise any deleterious impacts of inbreeding 
during creation of the transgenic strains. Our results revealed no 
significant difference between transgenic and wild-type mosquitoes 
for four of the fitness parameters measured, namely larval 
mortality, pupal sex ratio, hatch rate and longevity of blood-fed 
females. Some differences were identified for age at pupation, 
body size (wing length), fecundity and longevity of sugar-fed adults. 
However, these differences did not suggest transgenic inferiority. 
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Figure 4. Fecundity and fertility. Data represent comparisons 
between the control (KIL) and hemizygous transgenic strains (EV3P and 
EV3M). All data are pooled from 3 independent experiments. Different 
lower case letters (a, b, c) identify significant differences. (A) Wing 
length of adult females, reared for fecundity and longevity experiments 
with blood-meals. Boxplots show the interquartile range (open boxes), 
the median (horizontal line) and the maximum and minimum values 
(vertical lines). EV3M wing length is significantly longer than both KIL 
(p = 0.0003) and EV3P (p = 0.01 06). (B) Egg batch sizes, collected from 
single females showing the mean (open bars) and standard errors 
(vertical lines). The control KIL females produced significantly smaller 
egg batches than the transgenic strains EV3P (p = 0.0411) and EV3M 
(p<0.0001). EV3M egg batches were also significantly larger than EV3P 
(p = 0.0025). (C) Percent hatch rate from egg batches collected in (B). 
Boxplots show the interquartile range (open boxes), the median 
(horizontal line) and the maximum and minimum values (vertical lines). 
There are no significant differences in hatch rate between the control 
and either transgenic strain (p = 0.773). 
doi:10.1371/journal.pone.0088625.g004 

Whilst larval mortality did not differ between strains, the overall 
mean mortality of around 6% was significandy lower than that 
reported for An. gambiae in natural situations, where rates can be as 
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Figure 5. Adult longevity survival curves following mating and blood-feeding. Graphs represent comparisons between the control KIL 
(black) and hemizygous transgenic strains EV3P (green) and EV3M (red) and show Kaplan-Meier adult survival curves established with 3 day-old 
mated females. All survival curves represent pooled data from three independent experiments, each with 4 replicates of 25 adults for each strain. (A) 
Female survival in the absence of blood-feeding differs significantly between strains (p = 0.0004), with EV3P surviving significantly longer than either 
KIL (p = 0.002) or EV3M (p = 0.0002). (B) Female survival following blood-meals at days 0 and 5 of the experiment does not differ significantly between 
strains (p = 0.154). Significance levels are displayed on each graph (***, p<0.0001; *, ns; not significant). 
doi:1 0.1 371 /journal.pone.0088625.g005 



high as 92-98%. This can be due to varying factors such as the 
presence or absence of parasites, pathogens and predators within 
larval pools, inter- and intra-specific competition and abiotic 
factors such as temperature and rainfall [30,31]. With the 
exception of intra-specific competition, none of these macro- 
environmental variables are present in the laboratory setting. 
However, it has been reported that higher larval mortality rates 
from intra-specific or intra-strain competition may result in larger, 
more fecund adult females, resulting in no net fitness disadvantage 
[13]. 

No distortion in pupal sex ratios was found between strains since 
there were no significant departures from the expected 1 : 1 ratio of 
males to females. This is consistent with the facts that (a) the 
transgene is autosomal, (b) effector gene expression is predomi- 
nantly seen in blood-fed adult females and (c) fluorescence markers 
are regulated by the 3 xP3 promoter which is equally expressed in 
both sexes. It has been reported that carboxypeptidase mRNA is 
detectable in the pupal gut [32]. We therefore infer that mRNA 
and translated peptides from effector genes such as Vida3, which 
are expressed from the carboxypeptidase propter, are presumably 
also detectable in pupal tissues. However, our data show clearly 
that expression of the Vida3 effector peptide has no influence on 
pupal sex ratio. 

Age at pupation was remarkably consistent, with no significant 
differences between females of any strain and none between males 
of the KIL and EV3P strains. However, EV3M males were 
significantly slower to pupate than either KIL or EV3P males. 
Mating in An. gambiae occurs in male swarms, which await and 
then compete for, arriving females. Such swarms apparently lack a 
female-choice component, placing more emphasis on male 
competition [33]. It has also been suggested that slower 
development of transgenic mosquito larvae can delay the onset 
of male sexual maturity [12]. Thus, if EV3M males take longer to 
reach pupation and sexual maturity than their wild-type counter- 
parts, this could indicate a fitness disadvantage that might impact 
on transgene transfer to subsequent generations under field 
conditions. Conversely, EV3P males pupated at the same age 
and reached eclosion within a similar timeframe to KIL. 

There were no significant differences in wing length between 
strains for either males or females reared for the initial 



experiments. However, EV3M females subsequently reared for 
the assessment of fecundity, fertility and longevity following blood- 
feeding had significandy longer wing lengths. Wing length has 
been documented as a good indicator of body size for anopheline 
mosquitoes [34]. Furthermore, body size provides a good indicator 
of the blood-meal volume taken by females, which is a primary 
determinant of fecundity [35,36]. In one field study, smaller 
female mosquitoes, hand caught in Tanzania, were more often 
pre-gravid than their larger counterparts, requiring two or three 
blood-meals before completing the first gonotrophic cycle [35]. 
This study defined smaller females as having wing lengths of less 
than 3 mm whereas the laboratory strain females we describe here 
all had wing lengths in excess of 3 mm, indicative of optimized 
rearing conditions. Other studies of adult body size have reported 
that male An. gambiae preferentially select larger females, even 
when large numbers of males were competing for mates, thus 
suggesting that larger females are preferred due to their increased 
egg-carrying capacity [37]. Body size has also been reported to 
impact on male competitiveness during sexual selection in 
mosquito swarms [38]. This study suggested that the most 
successful males were of intermediate size (six times more likely 
to achieve copulation during one swarming event) although, within 
size groups, the mated males were larger than those which did not 
mate. It is therefore clear that body size of both males and females 
is an important determinant of reproductive success in An. gambiae. 
We therefore concentrate on the comparisons between EV3P and 
KIL females in the discussion of experiments involving blood- 
feeding below, where there were no significant differences in wing 
length. Since the use of transgenic mosquitoes in a disease control 
program would involve the release of homozygous males, it is 
important that such males are not compromised in body size and 
we have found no such evidence in studies of EVida3 homozygotes 
(unpublished data). 

Egg batches were significantly larger in transgenic line EV3M, 
which is entirely consistent with the longer wing length observed 
for these experiments. However EV3P females also produced 
larger egg batches than the wild-type control. Since there was no 
significant difference in wing length between EV3P and KIL 
females, we conclude that this transgenic strain is actually more 
fecund than the control. Regardless of egg batch sizes, hatch rates 
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were not significantly different between strains. Primigravida 
EV3P females could therefore be expected to produce more 
offspring than their KIL counterparts. 

Longevity experiments on mated females that went through two 
gonotrophic cycles were designed to determine whether upregula- 
tion of the transgene following naive blood-meals would impose a 
fitness cost. The An. gambiae carboxypeptidase promoter has 
previously been used to drive blood-meal inducible transgene 
expression in several An. stephensi transgenic lines. A fitness load 
was reported in two of these following naive blood-meals, both 
directly linked to the action of the transgenes themselves [6, 1 5] . 
For two different transgenic lines no significant effect on lifespan 
was observed following midgut transgene expression [15,24]. We 
observed no significant difference in longevity between transgenic 
and non-transgenic strains offered blood-meals. Thus, for EV3P 
females, where no significant differences to KIL were observed for 
wing length, hatch rate of progeny or lifespan, but for which egg 
batch sizes were significantly larger, there is the possibility of 
producing more progeny in a lifetime. We conclude that in the 
hemizygous transgenic mosquitoes tested here, neither the 
presence or the expression of the effector gene or fluorescence 
markers have a significant impact on larval mortality, pupal 
numbers and sex ratios, egg hatch rate or longevity following 
blood-feeding, and may in fact offer a fitness advantage since 
fecundity was increased. 

The greatest differences between strains were observed when 
assessing the longevity of non-mated and mated, non-fed adults. In 
the non-mated experiment, there was no significant difference 
between KIL and EV3M females but KIL males were significandy 
longer lived than EV3M males. Conversely, both EV3P males and 
females were significandy longer lived than their KIL counter- 
parts. This increase in lifespan, whereby EV3P females were 
significandy longer lived than KIL, was also observed for females 
that were mated and sugar starved. Thus, we identified significant 
differences between strains for adult longevity but, perhaps 
unexpectedly, the wild-type was not consistendy fittest. This 
suggests that there can be no a priori assumption of transgenic 
inferiority. We also observed some significant differences in 
longevity between reciprocal transgenic strains, whereby EV3P 
tended to be longer lived than their EV3M counterparts. This was 
unexpected since the two strains differ only in the route of 
inheritance of the transgene, from either the paternal or maternal 
parent. A similar situation was also evident from the data on age at 
pupation where male EV3M pupation was significandy slower 
than that of EV3P males. We believe this to be the first report of 
fitness disturbances associated with the route of inheritance of a 
transgene, although the mechanism remains unclear. It seems 
unlikely that such differences are a function of inbreeding, since 
both transgenic strains would have identical inbreeding coeffi- 
cients. We can therefore speculate that the route of inheritance of 
the transgene, rather than the transgene per se, can have pleiotropic 
effects on genes that influence fitness traits. This could arise from 
interactions with genetic determinants such as mRNA or proteins 
derived from the oocyte, or perhaps as a result of epigenetic 
mechanisms such as parental imprinting. It does, however, mean 
that the route of transgene inheritance must be considered in any 
biological or genetic characterization of transgenic strains destined 
for field release. 

The consequences of longevity variations in the field are likely to 
be complex. A longer lived female may undergo more gonotrophic 
cycles and thus increase reproductive fitness. Conversely, longevity 
may be less important to males, who can mate rapidly with 
multiple females even if lifespan is curtailed. There may even be a 
competitive advantage for those males that can mate early in life, 



although there is evidence to suggest a limit on mating capacity 
due to resource costs [39]. However, longer lived males will have 
greater opportunity to replenish resources and continue mating. If 
males are compromised in terms of reproductive fitness, perhaps 
because they take longer to pupate, and are less competitive early 
in life, then a reduced longevity would further compromise fitness. 
In addition, long-term survival in males is associated with high 
lipid content and larger body size [40]. However, a lower lipid 
content observed in laboratory-bred males (compared to natural 
males) could indicate a reduced reproductive potential [41]. 
Ng'habi et al. [38] also showed that larger mosquitoes were 
consistently longer-lived than intermediate sized males. They 
suggested that, although longer lived males have more opportunity 
to mate over their lifetime, this may not offset the fitness 
disadvantage their size confers since intermediate-sized males 
were found to be six times more likely to achieve copulation. 

The laboratory based fitness studies carried out here can not be 
truly representative of a field situation where mosquito larval 
stages will encounter increased resource competition, adults will 
forage for food (sugar and blood-meals) and both will be subjected 
to additional stress factors such as predation and pathogens. In all 
but our final longevity study (Figure 5B), females did not receive a 
blood-meal and so did not upregulate expression of the Vida3 
peptide or enter a gonotrophic cycle, all of which may affect 
longevity. Taking a blood-meal in the field risks exposure to 
pathogens, such as Plasmodium spp., which can upregulate immune 
responses, such as melanisation, and possibly decrease life 
expectancy [42]. The melanisation response is rarely seen in co- 
evolved Plamodium-znopheline combinations [43] and may 
therefore lead to a greater difference than expected between 
transgenic strains such as EVida3 and wild-type competitors. For 
transgenic refractory females, a reduction in average longevity due 
to upregulated immune responses [42] could be offset by reduced 
efficiency as a disease vector, offering an overall benefit. However, 
transgenic mosquitoes not previously exposed to pathogens or 
parasites in the laboratory, may elicit a greater immune response 
compared to native wild-type mosquitoes. Interestingly, transgenic 
malaria refractory lines were reported to be longer lived and more 
fecund when fed routinely on P. berghei-iniected blood in the 
laboratory [22]. We also noted that the virgin males used in our 
longevity studies lived longer than expected from general 
observation of stock populations within the same laboratory. This 
is consistent with the reports of Dao et al. [39] who showed that An. 
gambiae males exposed to females for six days had reduced 
longevity compared to virgin males. Fitness variations resulting 
from mating behaviour are likely to be even more pronounced in 
field situations where swarming is presumably more vigorous than 
in confined insectary conditions. 

Overall, the fitness parameters investigated in this study showed 
no evidence that the presence or expression of the Vida3 transgene 
or the associated fluorescence markers had a negative impact on 
fitness. Where significant differences between wild-type and 
transgenic strains were identified, they were not consistent across 
reciprocal hemizygous strains or between sexes. In addition, for 
some analyses, the transgenic lines outperformed their wild-type 
counterparts. Moreover, we identified variations in relative fitness 
associated with the route of transmission of the transgene (either 
from the male or female parent) that may need to be considered in 
the design and deployment of transgenic mosquitoes for field 
release. Such laboratory assessments of fitness are clearly necessary 
for initial characterization of newly generated transgenic strains. 
On their own, however, they would not be sufficient and would 
need to be supported by semi-field or field assessments of fitness 
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prior to any consideration for use in releases aimed at disease 
control. 

Materials and Methods 

Generation of Test Strain Mosquitoes 

Mosquitoes were maintained at 26°C ± 1°G and 80% RH in a 
1 2-hour light: 1 2-hour dark photoperiod. Trays (35 cm x25 cm) of 
wild-type KIT and homozygous transgenic EVida3 larvae (Li 
stage) were set up with equal numbers and reared following a 
standardised protocol until pupation. Pupae from both lines were 
sexed and equal numbers of the resulting adults put into each of 
three crosses to produce one KIT control population and two 
reciprocal hemizygous EVida3 strains (EV3P: EVida3 males with 
KIT females and EV3M: EVida3 females with KIL males). This 
protocol for generating experimental mosquitoes ensured the 
resulting wild-type and hemizygous strain individuals would be 
comparable. Adults were maintained on water and 10% glucose ad 
libitum. 4-7 day old adult mosquitoes were sugar starved for 18-24 
hours and females blood-fed on defribrinated horse blood (TCS 
Biosciences Ltd, Buckingham, UK) using a Hemotek membrane 
feeding system (Discovery Workshops, Accrington, UK) set at 
37°C. Polystyrene oviposition pots (5 cm diameter containing 
1 0 ml distilled H 2 0) were placed into each of the three cages 2 
days post blood-meal. Hatched larvae from crosses set up in this 
way were used in the experiments detailed below. 

Larval Survival, Age at Pupation and Sex Ratio 

Exactly 30 Li stage larvae were counted into each of 15 tubs 
(18.5 cmxlO cmx4.5 cm high) with 150 ml of distilled water 
containing Liquifry (Interpret Ltd., Dorking, U.K.). Tubs were fed 
daily on ground TetraMin fish flake (Melle, Germany) using a 
scaled down version of the standardised protocol below to ensure 
uniformity and water levels were kept constant throughout the 
experiment by addition of ddH 2 0 when necessary. Larvae were 
counted daily and dead larvae and pupae were removed. The 
number of male and female pupae per tub for each strain was 
recorded daily. Three independent replicates, derived from 
independendy reared adults were undertaken. 

Adult Longevity 

Experimental mosquitoes were generated using a standardised 
rearing protocol which ensured that both wild-type and hemizy- 
gous strain adults would be comparable. For wild-type KIL and 
both reciprocal hemizygous transgenic strains, L[ stage larvae 
were counted into trays at the standard density of 0.2 larvae per ml 
and distilled water containing 4 drops per litre Liquifry (Interpret 
Ltd., Dorking, U.K.). Thereafter a strict larval feeding regime was 
applied to all trays, using ground TetraMin fish flake or TetraMin 
baby (Melle, Germany) measured in a stainless steel microspoon, 
supplemented with Tetra Pond Sticks (Melle, Germany) at 
weekends. Trays of 200 larvae were fed l A flat spoon on day 2, 
1 flat spoon on day 4 and 2 flat spoons plus 1 2 pond sticks on day 
5. Water was topped up to the original level on day 5. From day 7 
onwards trays were given 1 flat spoon twice a day until all had 
pupated. Larval trays were rotated daily within the insectary to 
encourage synchronous growth and pupation. For initial longevity 
studies the resulting pupae were picked and sexed daily. Once 
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emerged, groups often adults of the same age and sex were placed 
in clear plastic pots (10 cm height x 11 cm diameter) with 10 
replicates of each sex per strain. For longevity studies on mated 
females, pupae of both sexes were collected and newly eclosed 
adults collected daily. Two days post-eclosion mosquitoes were 
sugar starved overnight. The following day, females of each strain 
were placed into eight separate 15 cm cages, 25 in 4 cages for 
non-fed samples and 30 in the remaining 4 cages for blood- 
feeding. Following the blood-meal any non-fed, partially-fed or 
excess females were removed to leave 25 in all cages. Fed cages 
were offered oviposition sites 2 and 3 days post-blood-meal and all 
cages were sugar starved for 3 hrs prior to offering a second blood- 
meal 8 days post-eclosion. Pots and cages were placed in the 
insectary in such a way as to minimise placement effects and 
rotated daily. Adults were maintained on water and sugar cubes 
available ad libitum. Dead mosquitoes were removed daily and 
counted for each population within each strain. The experiment 
ended when all mosquitoes in all populations were dead. Three 
independent replicates were undertaken using independendy 
reared adults for both experiments. 

Wing Lengths 

Adult mosquitoes from the wild-type strain and both reciprocal 
hemizygous strains were removed daily from pots and cages 
involved in the longevity studies and stored in 70% ethanol. These 
mosquitoes, reared following a standardised protocol (details 
above) had wing lengths measured as an indicator of body size. 
One wing was removed from each mosquito and measured from 
the distal end of the allula to the tip, excluding the fringe, using a 
dissecting microscope and an eye piece graticule calibrated using a 
1 mm stage micrometer ruler. Only intact wings were measured. 

Fecundity and Fertility 

Excess females reared for the second longevity study, which had 
been blood-fed 3 days post-eclosion, were used for fecundity 
studies. Fully engorged females were maintained on water soaked 
cotton wool and sugar cubes. At 3 days post-blood-meal single 
females were placed into 75 mmx25 mm (diameter), capped, soda 
lime glass specimen tubes (Fisher Scientific, Loughborough, UK), 
lined to 3 cm with filter paper and filled to 1.5 cm with distilled 
water. Tubes were examined daily. Once egg batches had been 
laid females were removed. Post-hatching, filter papers were slid 
out of the tube to collect eggs and stored in Petri dishes for 
counting. Total eggs laid per female (including any not collected 
on filter paper) were counted using a dissecting microscope, 
together with the number from each batch that had hatched. 
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